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Abstract 

We study the Higgs self couplings of a supersymmetric Eq model that has two 
Higgs doublets and two Higgs singlets. The lightest scalar Higgs boson in the model 
may be heavier than 112 GeV, at the one- loop level, where the negative results for 
the Higgs search at the LEP2 experiments are taken into account. The contributions 
from the top and scalar top quark loops are included in the radiative corrections 
to the one-loop mass of the lightest scalar Higgs boson, in the effective potential 
approximation. The effect of the Higgs self couplings may be observed in the pro- 
duction of the lightest scalar Higgs bosons in e + e~ collisions at the International 
Linear Collider (ILC) via double Higgs-strahlung process. For the center of mass 
energy of 500 GeV with the integrated luminosity of 500 fb _1 and the efficiency of 
20 %, we expect that at least 5 events of the lightest scalar Higgs boson may be 
produced at the ILC via double Higgs-strahlung process. 
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I. Introduction 



It is well accepted that one of the primary reasons to introduce the supersymmetry is 
is to elucidate the hierarchy problem of the Standard Model (SM) [1-3]. In the SM, 
the existence of the neutral Higgs boson is generally explained within the context of the 
origin of the electroweak symmetry breaking. The breaking of the electroweak gauge 
symmetry is induced by the self-interactions of the SM Higgs field, via the nontrivial 
vacuum expectation value. Without supersymmetry, the mass of the SM Higgs boson 
would acquire quadratically divergent loop corrections, requiring a cut-off scale at very 
high energy to confine the Higgs mass at the electroweak scale. 

If the supersymmetry is a symmetry of nature, the loop contributions of the parti- 
cles may be canceled by those of superpartners. However, since experiments have never 
observed yet any superparticles that have the same mass as quark or lepton, the su- 
persymmetry cannot be exact. In the presence of inexact supersymmetry, the one- loop 
corrections to the mass of the Higgs boson is not so small but as large as 30 % of the 
tree-level values, due to incomplete cancelation. 

There are a number of phenomenologically interesting supersymmetric models that 
embrace the SM. A characteristic of these supersymmetric standard models is the require- 
ment of at least two Higgs doublets, in order to give masses separately to the up-quark 
sector and the down-quark sector. The key role of these Higgs fields is essentially the same 
as the SM Higgs field: breaking of the electroweak gauge symmetry via self interactions. 

The simplest one is the minimal supersymmetric standard model (MSSM) [4], which 
has just two Higgs doublets and has been extensively studied. The Higgs sector of the 
MSSM can be extended by introducing Higgs singlet fields [5,6]. In the next-to-minimal 
supersymmetric standard model (NMSSM), there is one additional Higgs singlet. The 
motivation of introducing additional Higgs singlets is mainly to solve the fi problem of 
the MSSM [7]. In general, in terms of the vacuum expectation value of a neutral Higgs 
singlet, the \x parameter of the MSSM might be generated dynamically. 

Besides the MSSM and beyond, there are string-motivated supersymmetric models 
such as the supersymmetric Eq model, where the gauge symmetry is decomposed from 
Eq to 577(2) x U(l) x Ui(l) x U 2 (l). The two extra U(l) symmetries would be mixed 
to yield new linearly orthogonal combinations, U'(l) and U"(l). Thus, it is a rank-6 
supersymmetric model, with two extra neutral gauge bosons. Its Higgs sector has two 
Higgs doublets and two Higgs singlets. Phenomenology of the Higgs sector of this model 
has been investigated by several authors [8,9,10,11]. Some of us have elsewhere studied the 
possibility of detecting its neutral Higgs bosons at the Large Hadron Collider, assuming 
explicit CP violation in its Higgs sector [11]. 

The self interactions of the Higgs fields are induced by cubic and quartic couplings in 
the Higgs potential. Thus, in order to examine the Higgs phenomenology of a model, it 
is important to know the nature of the Higgs self interactions and study the possibility 
of measuring them. This kind of works have been performed some years ago in the SM 
and in the MSSM by Djouadi et al. [12,13]. 

In this article, we study the Higgs self couplings of the supersymmetric E 6 model that 
has been investigated in Refs. [8,9] which we will call the SUSY Eq model hereafter. We 
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study the possibility of discovering neutral scalar Higgs bosons of the SUSY E 6 model 
at the International Linear Collider (ILC). We calculate the cross section of the neutral 
Higgs bosons produced from double Higgs-strahlung process in e + e~ collisions. We find 
that, for reasonable parameter values, the cubic Higgs coupling of the lightest scalar 
Higgs boson is obtained to be smaller than the corresponding coupling of the SM Higgs 
boson. Thus, the production cross section of the lightest scalar Higgs boson via double 
Higgs-strahlung process is smaller than that of the corresponding SM Higgs boson. The 
radiative corrections to the mass of the lightest scalar Higgs boson due to the top and stop 
quark loops, combined with the negative experimental result for the neutral scalar boson 
search at LEP2 [14], suggest that the lightest scalar Higgs boson of the present model is 
heavier than 112 GeV. At the International Linear Collider (ILC) with the center of mass 
energy of 500 GeV, with the integrated luminosity of 500 fb _1 and the efficiency of 20 %, 
we expect that at least 5 events of the double Higgs-strahlung process would be observed. 



II. Higgs Potential 

The SUSY E 6 model has two Higgs-doublet superfields, Hi and H 2 , and two Higgs singlet 
superfields, A/l and M 2 . The superpotential of the model is given by [8.9] 

W « h t Q T H 2 t c R + XH 1 H 2 M l , (1) 

where t c R is the right-handed top-quark superfield, Q is the left-handed quark doublet 
superfield of the third generation, h t is the dimensionless Yukawa coupling coefficient for 
top quark, and A is a dimensionless coupling constant. The bottom-quark superfield is 
absent as we take only the top-quark superfield. Note that Af 2 is absent in the super- 
potential because the gauge symmetry of the model prohibits its coupling to the Higgs 
doublet superfields. 



The Higgs fields of the model are two Higgs doublets, Hi = (H®, H ) and H; 
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(H + , H®), and two Higgs singlets, Aq and N 2 . The fi parameter of the MSSM is dynam- 
ically generated in this model in terms of the vacuum expectation value of Aq. In terms 
of these Higgs fields, the Higgs potential at the tree level can be written as [9] 

U = m^H^ 2 + m 2 2 \H 2 \ 2 + ml\N x \ 2 + m\\N 2 \ 2 -(XAH^Nx + Yi.c.) 

+ |# 1 | 2 |Aq| 2 +|# 2 | 2 |Aq| 2 ] 

9 " mi 2 - \h 2 \y 




\2\ 



'2 



+ 9_ ll 
72 

+ Vl5C e (\H 1 \ 2 -\N 1 \ 2 + \N 2 \ 2 )} 2 , (2) 



^[Seim 2 + A\H 2 \ 2 - 5m 2 - 5\N 2 \ 2 ) 



where (i = 1, 2, 3, 4) are the mass parameters, g[ and g'[ are respectively the gauge 
coupling coefficients of U {!)' and U(l)", respectively, A is the trilinear soft SUSY-breaking 
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parameter with mass dimension, Cg = cos 6, and Sg = sin#, with 6 being the mixing angle 
between U(l)' and U(l)". 

Among the parameters in the Higgs potential, m; may be expressed in terms of other 
parameters by using the minimum equations that are obtained as the first derivatives 
of the Higgs potential with respect to the four neutral Higgs fields. After spontaneous 
breakdown of the electroweak gauge symmetry, we have two additional gauge bosons Z 
and Z for the extra U(l)' and U(l)", respectively, and seven physical Higgs particles: 
four scalar Higgs bosons, one pseudoscalar Higgs bosons, and a pair of charged Higgs 
bosons. The vacuum expectation values of the neutral Higgs fields are given as v\ = (Hi), 

v 2 = (H 2 ), Xi = (Ni), and x 2 = (N 2 ). We introduce tan/5 = v 2 /vi, and v = \Jv\ + v\ = 
175 GeV. 



III. Mass spectra 

Now, let us calculate the masses of the relevant particles. 



1. Masses of neutral gauge bosons and scalar top quarks 

In the SUSY Eq model, there are three neutral gauge bosons, Z, Z', and Z" arising from 
U(l), U(l)', and U(l)", respectively. Their square masses are given as the eigenvalues 
of a symmetric 3x3 matrix M G , the squared matrix. The elements of M G are given 
explicitly as 
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(3) 



where m 2 z 



v 2 (g\ + g 2 )/2, g\ and g 2 are the gauge coupling coefficients of U(l) and 
SU(2), respectively, S 2 g = sin(26 l ) and C 2 g = cos(26 l ). 
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A 3 x 3 orthogonal matrix U G can be employed to diagonalize M G , in order to obtain 
three eigenvalues, namely, m%, m 2 z ,, and m 2 z „. We note that the smallest eigenvalue, m|, 
is different from the squared mass of the neutral gauge boson in the SM. However, in the 
decoupling limit of the 27(1) symmetry from the extra U(l)' and U(l)", mz becomes the 
gauge boson mass of the SM. We sort them such that mz < mz> < mz", where the mass 
of the lightest neutral gauge boson should be equal to the mass of the existing neutral 
gauge boson, 91.2 GeV. 

The mixings among them can be parameterized by three mixing angles, namely, a% 
between Z and Z', a 2 between Z' and Z", and «3 between Z and Z", respectively. These 
mixing angles may be expressed in terms of the elements of U G as 



fug' 

ati = arctan — 

\U 2 2. 



a 2 



arcsm 
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ug, 



(4) 



We note that there are strong experimental constraints on the masses of the extra 
neutral gauge bosons and on the mixing angles between the extra gauge bosons and 
the SM Z boson. The results of direct search in pp collisions at the Fermilab Tevatron 
constrain that rnz', rriz>> > 800 GeV, and the electroweak precision measurements at 
LEP2 constrain that \a±\, \as\ < 3 x 10 -3 [15]. Since the mixings between 27(1) and 27(1)' 
or between U(l) and 27(1)" are very small, the extra gauge symmetries are more or less 
decoupled from SU(2) x 27(1). 

Next, we calculate the scalar top quark masses at the tree level. They are obtained as 



- (m 2 Q + n4) + mf + -^rn 2 z c< 
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where rriQ and ttit are the soft SUSY breaking masses, 



(7) 



v 2 g 2 /2, m t = h t V2 is the 



top quark mass, A t is the trilinear SUSY breaking parameter with mass dimension, G t 
is the D-term contribution from the extra 27(1)' and 27(1)", and X t is the mixing term 
between t\ and t%. 
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2. Masses of pseudoscalar and scalar Higgs bosons 

Next, we evaluate the one-loop corrections to the mass spectra. The corrections are 
calculated by inserting the relevant tree-level masses into the one-loop effective potential 
a la Coleman and Weinberg [16]. The top quark and scalar top contributions to the 
one-loop effective potential is given by 
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niMf 
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log 
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where A is the renormalization scale of the modified minimal substraction scheme. The 
summation is over I = ti, t 2 , t and rii = 6 for the scalar top quarks and ri\ = —12 for the 
top quark. 

There are four neutral Higgs fields in the SUSY E§ model, arising from two Higgs 
doublets and two Higgs singlets. Among four complex components of these neutral Higgs 
fields, three of them are gauged to generate masses of Z, Z' ', and Z"; and the remaining 
one complex component becomes the physical pseudoscalar Higgs boson of the model. 
The mass of the pseudoscalar Higgs boson at the one-loop level is given as 
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where m 2 A0 is obtained from the tree-level Higgs potential Vo, m A i is obtained from the 
one-loop Higgs potential Vx, tana = v sin 2/3/(2xi), and 
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The second derivatives of the Higgs potential with respect to the real components of 
these four neutral Higgs fields give us the mass matrix for four scalar Higgs bosons. The 
4x4 symmetric squared mass matrix M for the scalar Higgs bosons can be decomposed 

as 



M = M° + M 1 



(12) 



where M° is obtained from the tree-level Higgs potential Vo, and M 1 from the one-loop 
Higgs potential V\. The explicit formulas for M°- (z, j = 1, 2, 3, 4) are obtained as 
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and the explicit formulae for Ml (2, j = 1,2,3,4) are obtained by inserting the Higgs field 
dependant masses of the top quark and scalar top quarks at the tree-level into V 1 as 
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The eigenstates of M are four physical scalar Higgs bosons of the SUSY Eq model, denoted 
as Si (i — 1, 2, 3, 4), and the corresponding eigenvalues are their squared masses, denoted 
as m|. . The masses of these four scalar Higgs bosons are ordered as m Sl < m S2 < m Si < 

For the numerical analysis, we establish the parameter space as follows: 1 < tan (3 < 
30, < A < 0.83, 10 < A < 400, < 9 < tt/2, 100 < x 1: x 2 < 1500 GeV, 100 < 
Trig, rriT, A t < 1000 GeV. Furthermore, we put m t = 175 GeV, and we assume that the 



8 



scalar top quarks are heavier than the top quark. We set the lower bound on the effective 
fi e = Xxi parameter to be 150 GeV, in order to take into account the experimental lower 
bound on the lighter chargino mass at the LEP2 experiments. 

We also take into account the experimental results for the Higgs search at the LEP2 
experiments. Recently, the collaborations of the LEP2 experiments have reported the 
model-independent upper bound on the Higgs coupling coefficient to a pair of Z bosons 
at 95 % confidence level [14]. The result may be interpreted for the coefficient of the SUSY 
Eq model. The coupling coefficient of Sj to a pair of Z bosons of the model, normalized 
by the corresponding SM coupling coefficient, can be written as 

G ZZSl ~ cos (30 u + sin (30 2i , (20) 

where are the r/'-th element of the 4x4 orthogonal matrix that diagonalizes the mass 
matrix for the neutral scalar Higgs bosons. 

For a given set of parameter values, we evaluate m s . and GzzSi- m this way, using 
Monte Carlo method, we explore 10 5 points of the established parameter space. In Fig. 
1, we show the result of evaluations of m Sl and Gzzs x - The SM coupling coefficient is 1.0 
since it is also normalized by itself. We find that GzzSi ls larger than 0.9 for the most 
of the parameter space, which tells us that S\ is more or less equivalent to the SM Higgs 
boson. Since GzzSi = 1; 2, 3, 4) satisfy a sum rule of J2i=i G\ zs . ~ 1, this implies that 
the pair of Z bosons couples nearly exclusively to S\ in most of the parameter space and 
the couplings to the other heavier scalar Higgs bosons are negligibly small. 

The result in Fig. 1 also shows that ms 1 is in the range of 112 to 142 GeV at the 
one- loop level. The upper bound on is determined primarily by the maximum value 
of the SUSY breaking scale and the LEP2 constraints on the SM Higgs coupling coefficient 
to a pair of S bosons. The masses of heavier scalar Higgs bosons are 135 GeV < ms 2 < 
897 GeV, 800 GeV < m s . A < 1155 GeV, and 1033 GeV < m Si < 2828 GeV. Also, the 
masses of the pseudoscalar Higgs boson of the SUSY E 6 model at the one-loop level are 
estimated to be 120 < rriA < 2828 for the parameter values we consider. 

The coupling coefficient Gzzs x is crucial to calculate the production cross section for 
the Higgs-strahlung process, <r(e + e~ — ► Z — > ZS\). In Fig. 2, we show the production 
cross section of Si via Higgs-strahlung process in e + e~ collisions with s/s = 500 GeV, 
where the corresponding SM cross section is also shown in a solid curve, as a function of 
the SM Higgs boson mass. The values of the relevant parameters are the same as in Fig. 
1, chosen randomly in the established parameter space. 

One may notice in Fig. 2 that most of the 10 5 points are distributed quite close to 
the solid curve. In other words, Si of the SUSY Eq model behaves almost the same as 
the SM Higgs boson in the Higgs-strahlung process in e + e~ collisions. One may also 
notice in Fig. 2 that the absolute minimum of the production cross section is about 19 
fb. That is, parameter values, the production cross section is larger than 19 fb for any. 
It is worthwhile to compare this absolute lower bound on the production cross section 
with the corresponding results of other models. We have studied the same subject in the 
next-to MSSM (NMSSM) [17]. We have calculated that the absolute lower bound on the 
production cross section of the lightest scalar Higgs boson via Higgs-strahlung process in 
the NMSSM is about 15 fb (see the dashed curve of Fig. 3a in Ref. [17]). Thus, the 
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absolute lower bound on the production cross section of the SUSY Eq model is in fact 
larger than the corresponding value in the NMSSM. 

It is usually anticipated that, for a supersymmetric model with several Higgs singlets, 
the production cross section for the lightest scalar Higgs boson tends to decrease as the 
number of Higgs singlets increases, because the probability of production should be shared 
with other heavier scalar Higgs bosons. In this respect, the lower bound on the production 
cross section of the SUSY Eq model should be smaller than that of the NMSSM, since 
the SUSY E 6 model has one more Higgs singlet than the NMSSM. However, the result is 
opposite to the usual anticipation. We think that the main reason for this result is that 
Si of the SUSY Eq model is more similar to the SM Higgs boson than the lightest scalar 
Higgs boson of the NMSSM, at least in the Higgs-strahlung process in e + e~ collisions. 

IV. Production of scalar Higgs bosons via double 
Higgs-strahlung process 

The cubic and quartic couplings of Higgs bosons contribute when a process involves mul- 
tiple Higgs bosons. Typically, the double Higgs production process in e + e~ collisions 
depends on the cubic couplings of Higgs bosons and the quartic couplings of Higgs-Higgs- 
gauge-gauge bosons. For practical simplicity, we only consider double S\ productions: 
e + e~ — > ZSiS\. The double productions of heavier scalar Higgs bosons are kinematically 
suppressed and therefore the production cross sections are very small. The relevant Feyn- 
man diagrams for the double S\ production are shown in Fig. 3. Notice that the quartic 
coupling and the cubic coupling are present in Fig. 3. 

In principle, all of the four scalar Higgs bosons of the SUSY Eq model may participate 
in Fig. 3 as the intermediate particle. The contributions of the heavier scalar Higgs 
particles as the intermediate particle depend on two coupling coefficients: GzzSi a t the 
one end and Gs i s 1 Si a ^ the other end (i = 2, 3, 4). However, as we have observed before, 
GzzSi = 2, 3, 4), the coupling coefficients of heavier scalar Higgs bosons to the Z boson 
pair, are negligibly smaller than GzzSi- Thus, the contributions of the heavier scalar 
Higgs particles as the intermediate particle may be neglected, and we only consider Si 
as the intermediate particle in Fig. 3 in our numerical analysis. Consequently, only Si 
participates in Fig. 3. 

In order to calculate the production cross section of the lightest scalar Higgs boson, 
we have to know its cubic and quartic coupling coefficients in the SUSY Eq model. The 
cubic coupling coefficient of Si can be written as 

Gs^Si = GsiSiSi + G s t SiSi ( 21 ) 

where G^ , s . s . is the tree-level coefficient and Gg. s . s . is the one-loop correction. Explicitly, 
the tree- level cubic coupling coefficient of Si is given as 

G% SiSi = 12A 2 [(0* + 2 2l )0 3l xi + OiAOl + 3 >cos/3 + 2i (0 2 u + 3 2 >sin/3] 

+Mgi + gDiol - o 2 ^{o u cos p - o 2l sin/3) - 6m ^ oSm2a OAo 3t 
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where Go <? c the D-term contribution of the U(l)' and Go a <? the D-term contribution 
of the £7(1)" at the tree level. 

The one-loop correction to the cubic coupling of the scalar Higgs boson is given as 
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The corresponding SM cubic coupling coefficient, Ghhh, may be obtained from the SM 
Higgs potential /x 2 0V/2 + A(<^ <fr) 2 / '4 as Qml/v, where v = ((f)) and is the mass of SM 
Higgs boson. Conventionally, the SM cubic coupling coefficient is normalized such that 
[13] Ghhh = Ghhh ■ u /(2w|). In Fig. 3, we also plot G\ hh as a function of rxih (the solid 
curve). 

In order to compare the cubic coupling coefficient of the SUSY E$ model with the SM 
cubic coupling coefficient, we may normalize Gs^Si m ^ ne un ^ °f 2m 2 z /v as [v/(2m z )]Gs i s i s i - 
In Fig. 4 we plot the square of the cubic coupling coefficient, normalized by 2m 2 z /v, against 
msx, by varying the values of the relevant parameters within the established space. We 
find that the cubic coupling coefficient of the SUSY E§ model is smaller than the corre- 
sponding SM coefficient. This behavior of the SUSY Eq model is similar to the MSSM, 
where the cubic coupling coefficient of the lightest scalar Higgs boson is smaller than the 
corresponding coefficient of the SM Higgs boson [18]. On the other hand, it is different 
to the two Higgs doublet model, where the SM cubic coupling coefficient may be smaller 
than the cubic coupling coefficient of the lightest scalar Higgs boson [19]. We also find 
that the cubic coupling coefficients of other scalar Higgs bosons, Gs^St {i — 2,3,4) are 
smaller than the SM cubic coupling coefficient. 

Next, let us consider the quartic coupling coefficient of Higgs- Higgs-gauge-gauge bosons 
of the SUSY E 6 model. When normalized by the corresponding SM quartic coupling co- 
efficient, GzzSiSi °f Si pair to Z boson pair is given as 



G 



Z Z S{Si 



o 2 



o 2 

w 2i 1 



(26) 



if the mixing between Z boson and extra neutral gauge bosons is neglected. Since the 
mixing angles \a 2 \ and | CK3 1 are very small, the neglection is practically justified. 

Now, we are ready to calculate the differential cross section for the double Higgs- 
strahlung process in e + e~ collisions of the lightest scalar Higgs boson of the SUSY E 6 
model. It is obtained as 



da 



dx\dx<2 



e e 



Z -> ZSxSx 



V2G 3 f 



v 2 + a 2 



384vr 3 s 



Hz)'' 



Z 



(27) 



where Xj = 22^/ \fs is the scaled energy of the i-th Si, with E^ being the energy of the 
z-th Si and \fs being the cm. energy of e + e~ system (i = 1,2), v e = — 1 + 4 sin 2 9w and 
a e = — 1 are respectively the vector and axial- vector Z charges of the incoming electron, 
\iz = tti% /s is the square of the reduced Z boson mass, and Z is given by 
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ZiG 
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2 n 2 
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Z V2 G ZZS X G S^S^ 



Zi^G zzsfi ' s^SxSxG ' zzsxSx + 2 2 3G ZZSi Gzzs 1 s 1 , (2i 
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x\ — X2 is the scaled energy of the Z boson, y\ = 1 — x±, y% — 1 
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Hh = m 's 1 / s is the square of the reduced mass of Si, and \inz = I^h — Hz- 

Then, we integrate the differential cross section to obtain the total production cross 
section for the double Higgs-strahlung process of Si, a(ZS\S\) as a function of ms 1 . The 
values of the relevant parameters are tha same as in Fig. 1, 10 5 points in the parameter 
space. The result is shown in Fig. 6, where we take yfs = 500 GeV, the proposed cm. 
energy of the first stage of the ILC (ILC-500). 

For the sake of comparison, we also plot the corresponding SM cross section in Fig. 
5. One may notice that most of the points are distributed close to the solid curve. This 
implies that Si behaves very much alike the SM Higgs boson with respect to the double 
Higgs-strahlung process. Also, the absolute minimum of the production cross section of 
the lightest scalar Higgs boson via the double Higgs-strahlung process is about 0.05 fb at 
the ILC-500. In other words, the SUSY E 6 model predicts that a(ZSiSi) is larger than 
0.05 fb, whatever the parameter values. Therefore, we expect that the ILC-500 would 
produce at least 5 Si events for the SUSY E 6 model via the double Higgs-strahlung 
process, if the ILC-500 has the integrated luminosity of 500 fb -1 and the efficiency of 20 
%. 

V. Conclusions 

We have studied a supersymmetric Eq model, which has two Higgs doublets and two 
Higgs singlets. In particular, we have studied the production of the lightest scalar Higgs 
boson in e + e~ collisions, via Higgs-strahlung process and double Higgs-strahlung process, 
for a reasonably established parameter space, at the one-loop level by considering the 
contributions from top and scalar top quark loops. 

At the one-loop level, the mass of the lightest scalar Higgs boson of the SUSY E% model 
is estimated to be between 112 and 142 GeV. Thus, its mass is comparable to the SM 
Higgs boson. Not only its mass but also its production cross sections in e + e~ collisions via 
Higgs-strahlung process and double Higgs-strahlung process are found to be quite similar 
to those of the SM Higgs boson. For the most part, this similarity may be attributed to 
the experimental constraints on the large masses of the extra neutral gauge bosons and 
the tiny mixings between them and the Z boson, the electroweak precision measurement 
at the LEP2 experiment and the direct search in pp collisions at the Tevatron. Since the 
extra neutral gauge bosons are heavier than 800 GeV and their mixing angles with the Z 
boson are less than 3 x 10~ 3 , they are practically decoupled from Z boson, and hence from 
Si in the Higgs-strahlung process and double Higgs-strahlung process. In other words, 
we could assume without difficulty that only the Z boson is involved in these processes. 
We also find that the same experimental constraints set the lower bound on x\ and X2, 
the vacuum expectation values of two Higgs singlets, as about 1265 GeV. 

Our study on the Gzzs x coupling coefficient, the production cross section <r(e + e~ — > 
Z — > ZZSi) via Higgs-strahlung process, the Gs 1 s 1 Si cubic coupling coefficient, Gzzs x s x 
quartic coupling coefficient, and the production cross section <r(e + e" — > Z — > ZS\S\) via 
double Higgs-strahlung process show clearly the similarity of Si of our model with the SM 
Higgs boson. Thus, the contributions of heavier scalar Higgs bosons are nearly negligible 
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in these processes. Consequently, the experimental constraints on the extra neutral gauge 
bosons make the Higgs sector of the supersymmetric E 6 model with two Higgs doublets 
and two Higgs singlets very similar to the SM Higgs sector, as far as Si is concerned. 

The absolute lower bound on the cross section for Si production in e + e~ collisions via 
Higgs-strahlung process at the ILC-500 is about 19 fb, and the absolute lower bound on 
the cross section for Si Si production in e + e~ collisions via double Higgs-strahlung process 
at the ILC-500 is about 0.05 fb. Assuming the integrated luminosity of 500 fb -1 and the 
efficiency of 20 % for the ILC-500, at least 5 events of Si of our model might be explored 
at the ILC via the double Higgs-strahlung process. 
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110 115 120 125 130 135 140 145 

m s1 (GeV) 

FIG. 1: The squared coupling coefficient G 2 ZZSi against ms 1 , normalized by the SM coef- 
ficient, for 10 5 sets of parameter values that vary within the parameter space established 
as 1 < tan/5 < 30, < A < 0.83, 10 < A\ < 400, < 9 < vr/2, 100 < x u x 2 < 1500 
GeV, 100 < rriQ, rriT, A t < 1000 GeV. The squared SM coefficient, normalized by itself, 
is a solid line at 1.0. Notice that most of the points are distributed close to the solid line, 
and that m Sl is between 112 and 142 GeV. 
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FIG. 2: The production cross section a(ZZSi) against m^, via Higgs-strahlung process 
in e + e~ collisions with -^/i = 500 GeV, for the same 10 5 sets of parameter values as in 
Fig. 1. The solid curve is the corresponding SM cross section as a function of the SM 
Higgs boson mass. 
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Si 



FIG. 3: Feynman diagrams for double Higgs-strahlung process, e + e~ — > Z — > ZS\S\. 
Notice that the Higgs cubic coupling and the Higgs-Higgs-gauge-gauge quartic coupling 
are present. 
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(^S1S1Sl) 
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m s1 (GeV) 

FIG. 4: The squared coupling coefficient G^ lSlS against ms 1 , normalized by 2m 2 z /v, 
for same 10 5 sets of parameter values as in Fig. 1. The solid curve is the squared SM 
coefficient, normalized by 2m z /v , as a function of the SM Higgs boson mass. 
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FIG. 5: The production cross section a(ZSiSi) against m Sl , via double Higgs-strahlung 
process in e + e~ collisions with y/s = 500 GeV, for the same 10 5 sets of parameter values 
as in Fig. 1. The solid curve is the corresponding SM cross section as a function of the 
SM Higgs boson mass. 
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